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IV 


ABSTRACT 


The  shallow-water  environment  poses  many  obstacles  to  submerged  and  surface 
operations.  Not  least  of  them  are  the  obstacles  posed  to  acoustic  detection  of  quiet 
contacts.  The  presence  of  Internal  Solitary  Waves  (ISWs)  within  this  environmet  induce 
acoustic  variability  and  phase  fluctuations,  which  limit  signal  processing  capability  and 
further  complicate  operations  in  this  environment.  However,  a  better  understanding  of 
how  ISWs  affect  sound  propagation  may  lead  us  to  take  advantage  of  some  of  the  effects 
of  the  ISWs,  such  as  how  ISWs  focus  acoustic  energy  into  certain  modes  allowing 
detection  of  quiet  sources  above  background  noise.  An  understanding  of  the  phase 
fluctuation  may  also  lead  to  better  processing  algorithms.  While  multiple  numerical 
simulations  have  been  conducted,  which  allow  for  the  effects  of  ISWs  to  be  intuited  from 
them,  an  analytical  model  that  can  predict  the  effects  on  acoustic  propagation  by  the 
physical  parameters  of  the  ISWs  had  not  been  produced  until  Professor  John  A.  Colosi  of 
the  Naval  Postgraduate  School  developed  a  single  scattering  analytical  model  of  those 
effects. 

In  this  thesis,  the  Colosi  Single  Scattering  model  will  be  compared  to  an  accepted 
3-D  PE  model  in  order  to  further  validate  the  model.  In  the  comparison,  numerous 
simulations  from  both  models  have  been  used  to  determine  how  closely  the  Single 
Scattering  model  can  predict  the  vertical  mode  coupling  effects  of  ISWs. 
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I.  INTRODUCTION 


Shallow-water  Internal  Solitary  Waves  (ISW)  have  been  observed  in  a  variety  of 
loeations  ineluding:  the  Massaehusetts  Bay,  the  New  York  Bight,  Gulf  of  California, 
Andaman  Sea  offshore  Thailand,  the  Australian  North  West  Shelf,  the  Sulu  Sea  between 
the  Philippines  and  Borneo,  off  the  eoast  of  Portugal,  off  Hainan  Island  in  the  South 
China  Sea,  and  off  the  Strait  of  Gibraltar  in  the  Alboran  Sea  (Zhou,  Zhang,  and  Rogers, 
1991).  These  shallow-water  ISWs  have  soliton  like  eharaeterizations  sueh  as:  well- 
defined  wavelengths,  propagation  in  paekets  towards  the  shoreline,  higher  than  linear 
group  veloeity,  and  deerease  in  wavelength  and  amplitude  toward  the  rear  of  the  paeket 
(Zhou,  Zhang,  and  Rogers,  1991).  The  solitons  are  caused  by  tidal  driven  flows  over 
unusual  bathymetry,  such  as  sills  and  edges  of  continental  shelves,  and  occur  during  the 
summer  months  when  there  is  a  strong  thermocline  near  the  surface  and  the  layer  is 
highly  stratified  (Zhou,  Zhang,  and  Rogers,  1991).  Figure  I  shows  the  effects  of  a  soliton 
packet  on  the  sound  speed  profile  (ssp)  within  the  water  column.  These  perturbations  to 
the  ssp  can  cause  unusually  high  transmission  loss  due  to  vertical  mode  coupling.  This 
transmission  loss  has  a  resonance  like  behavior  that  is  a  function  of  acoustic  source 
frequency  and  effective  width  of  the  soliton,  where  the  effective  width  of  the  soliton 
refers  to  the  distance  traveled  by  the  acoustic  energy  passing  through  each  individual 
soliton  (Zhou,  Zhang,  and  Rogers,  1991;  Colosi,  2008). 
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Figure  1.  A  Three  Soliton  Packet  Supper  Imposed  on  the  SSP  of  a  100  m  Deep 

Sound  Channel 

Prior  to  Professor  John  Colosi  of  the  Naval  Postgraduate  School  developing  the 
Single  Scattering  Analytical  Model  for  2-D  acoustic  mode  propagation  through  a  soliton, 
there  were  no  analytical  descriptions  of  what  effects  solitons  have  on  vertical  mode 
coupling  (Colosi,  2008).  An  analytical  model,  which  can  mathematically  model  the 
effects  of  solitons  on  mode  coupling,  is  beneficial  to  improving  processing  algorithms. 
This  paper  seeks  to  validate  the  Single  Scattering  Analytical  Model  by  direct  comparison 
of  its  results  with  a  3-D  Parabolic  Equation  Model.  A  Parabolic  Equation  model  is  used 
as  a  benchmark  in  this  analysis,  since  this  type  of  model  has  been  proven  to  produce 
results  consistent  with  real  world  experiments  (Smith,  2001).  The  results  will  be 
compared  for  propagation  through  both  a  single  soliton  and  a  soliton  packet  for  multiple 
orientations  of  the  solitons.  The  comparison  will  be  based  on  the  ratio  of  modal  energies 
(RME).  The  RME  is  a  measure  of  energy  coupled  into  or  out  of  a  particular  mode.  The 
resonance  behavior  exhibited  by  the  mode  coupling  is  compared  as  well  by  predicting  the 
resonance  location  of  the  PE  model  using  the  analytical  model. 
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II.  HOW  THE  PE  MODEL  WORKS 


A.  INTRODUCTION 

In  this  chapter,  a  brief  description  of  the  3-D  implementation  of  the  Monterey 
Miami  Parabolie  Equation  (3DMMPE)  model  is  given  (Smith  and  Colosi,  2008).  The 
description  shall  cover  the  pertinent  material  from  the  SWAM99  article  Professor  Smith 
published  in  the  Journal  of  Computational  Acoustics  (JCA)  (Smith,  2001),  and  the 
extension  to  3-D.  The  first  part  of  the  discussion  covers  how  the  parabolic  equation  was 
developed.  Next,  a  discussion  of  the  split-step  Eourier  algorithm  is  given.  Easily,  the 
method  for  modal  deeomposition  is  explained. 

B,  DEVELOPMENT  OF  THE  PARABOLIC  EQUATION 

In  this  seetion,  the  development  of  the  3DMMPE  is  discussed,  as  well  as  aspeets 
that  are  unique  to  this  particular  PE  model.  This  seetion  closely  follows  the  SWAMM99 
article  for  the  development  of  the  PE  (Smith,  2001).  This  referenee  is  presented  to 
prevent  exeessive  duplieate  references  in  the  following  section.  Unlike  the  SWAMM99 
artiele,  the  PE  is  developed  in  terms  of  Cartesian  coordinates  instead  of  cylindrical 
coordinates. 

Eirst,  the  time-harmonie  aeoustic  field  is  defined  in  Cartesian  coordinates  as 

P{x,y,z,t)  =  p{x,y,z)e““‘ .  (1) 

Substitution  into  the  wave  equation  results  in  the  Helmholtz  equation, 

V V(U  T,  z)  +  kin  (x,  y,  z)p{x,  y,  z)  =  -A7tP^5{z  -  zj  (2) 

The  value  is  the  reference  sound  speed  and  is  defined  as  k^  =  col c^,  and  n  is  the  index 
of  refraction  defined  as  n  =  c^l  c{x,y,z) .  The  value  is  the  reference  acoustic  sound 

speed,  which  is  picked  so  that  all  values  of  are  near  unity.  The  acoustic  sound  speed 
c(x,y,z)  contains  all  of  the  features  of  the  environment  including  the  sound  speed 
perturbation  caused  by  the  solitons.  The  term  on  the  right  hand  side  of  equation  is  the 
source  term  which  approximates  the  souree  loeated  at  x  =  0,  y  =  0,  and  z  =  z^.  Setting  the 
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right  hand  side  to  zero  allows  the  homogenous  Helmholtz  equation  to  be  solved  for  the 
aeoustie  pressure.  Rewriting  the  Helmholtz  equation  in  terms  of  the  partial  derivatives 
yields 


p  d^p  d^p 
dx^  dy^  dz^ 


+  kln^ p  =  0  . 


(3) 


Introdueing  the  operators 

P  =—. 
&' 


1  5' 


1  5' 


u  = 


allows  the  Helmholtz  to  be  rewritten  as 


V  = 


V  +  U  +  £  + 


kl  dy^ 

l^p  =  0 . 


and  £  =  n 


(4) 

(5) 


Introdueing  the  operator  ={u+v  +  s  +  further  simplifies  the  Helmholtz  equation 


to 

(6) 

Equation  6  ean  be  faetored  to 

(K,  +  Ik'-QPliP., -ik'-Q,, )  P  +  ik.  [/>„,,  a,  ]  P  =  0 .  (7) 

where  it  represents  the  factorization  of  the  field  into  incoming  and  outgoing  waves 
(Smith,  2001;  Smith,  2005).  If  the  commutator  {\_Pgp,QopV)  is  taken  to  be  zero  for  layered 

media,  and  the  substitution  p  =  Q^p^^  is  made,  the  outgoing  wave  then  satisfies 

Pop^=iKQop^-  (8) 

Defining  'F  in  terms  of  the  PE  field  function  y/  such  that  'F  =  y/e'""^ ,  then  the  outgoing 
field  function  (Equation  8)  can  be  written  in  terms  of  the  evolution  of  the  PE  field 
function  in  the  x-direction,  i.e., 

^  =  ikoiQop  -  =  -iKHop¥  (9) 

The  acoustic  pressure  field  can  then  be  determined  from  the  PE  field  function  by 

p{x,  y ,  z)  =  P^  A—Qkp^y/{x,  y ,  z)e  .  (10) 

V  r 
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Here,  refers  to  the  source  acoustic  pressure  at  a  reference  range  ,  and  r  is  the 
horizontal  slant  range  r  =  +y^ 


C.  DESCRIPTION  OF  THE  SPLIT-STEP  FOURIER  ALGORITHM 


Now  that  an  expression  for  the  evolution  of  the  PE  field  function  in  range  (x- 
direction)  has  been  developed  in  Equation  9,  all  that  is  needed  is  an  algorithm  to  solve  for 
y/  .  Equation  9  can  be  rewritten  in  the  form 

^(r  +  Ar)  =  (h(r)^(r)  ,  (11) 

where  (h(r)  represents  the  propagator  of  y/{r) .  Eor  small  range  steps  Ar ,  the  propagator 
can  be  approximated  by  (Smith,  2005) 

=  (12) 

To  implement  the  SSE  algorithm  the  ,  and  therefore  the  ,  must  be 

separable.  Specifically,  the  u  and  v  operators  need  to  be  separated  from  the  s  operator. 
In  the  3DMMPE,  this  is  accomplished  with  the  Thompson  and  Chapman  (Smith,  2005) 
approximation  for  the  .  This  approximation  allows  the  to  be  represented  by 

Qop  —  Vl  +  M  +  V  +  £■  «  Vl  +  M  +  V  +  yj\  +  s  “  1  =  1  “  Tgp  “  Ugp  ,  (13) 


where 

T^p  =\-y/\-u+v  and  =\-^\  +  s  .  (14) 

Recall  from  Equation  9  that  =  ^-Qop,  which  when  substituted  into  Equation  12 
yields 

^  ^-iKAr(Q„p  -1)  ^  +C/„p  )  ^  ^  J 

Using  the  Baker-Campbell-Hausdorff  expansion  (Bellman,  1964),  Equation  15  can  be 
rewritten 


(16) 


Substituting  Equation  16  back  into  Equation  1 1  yields 


ik„—U„Jr+hr)  ^  ik^—U„Ar) 

^(r  +  Ar)  =  e  ^  e  ^  'e  ^  ■ 


(17) 
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The  is  a  scalar  in  z-space  and,  therefore,  is  a  diagonal  matrix  that  can  be  multiplied 

with  .  The  is  a  differential  operator  in  terms  of  z  and  y,  so  that  different  depth 

eigenfunctions  are  coupled.  As  a  result,  the  is  not  a  diagonal  operator  in  z-space  and 

cannot  just  be  multiplied  with  the  y/{r)  term.  The  operation  of  on  ^(r)  in  (y,z)- 

space  would  then  be  a  time  and  memory  intensive  process.  An  alternative  approach 
would  be  to  take  the  Fourier  Transform  of  y/{r)  and  then  operate  on  it  in  (k^,  k^j-space . 

The  benefit  of  performing  the  operation  in  (k^,A:^)-space  is  that  is  a  diagonal  matrix 

and,  therefore,  is  a  scalar  multiplication  operation  with  yjr{r) .  Here,  is  equal  to 

The  inverse  Fourier  Transform  of  y/{r)  can  then  be  taken  before  performing  the  scalar 
multiplication  operation  with  .  This  is  exactly  the  approach  taken  in  the  Split-Step 
Fourier  Algorithm,  shown  symbolically  below. 


kl+kl 


(18) 


y/{x  +  lxx,y,z)  =  e 


-K-u,,, 

o  2  "i 


(x+Ax) 


IFFTle 


-ik„&xT(k^,k) 


FFT 


o  ^  <>1 

o  Z 


(X) 


y/{x,y,z) 


(19) 


D,  METHOD  FOR  MODAL  DECOMPOSITION 

In  this  section,  the  method  used  to  extract  the  modal  amplitudes  ( ^„ )  from  the 

acoustic  pressure  field  computed  by  the  PE  model  is  discussed.  This  discussion  includes 
why  the  normal  mode  shape  can  be  used  as  the  basis  set  for  the  wave  equation.  Errors 
introduced  by  the  Thompson  Chapman  approximation  will  also  briefly  be  discussed. 

Smith  used  a  separation  of  variables  technique  to  obtain  the  depth  separable  wave 
equation  (Smith,  2005).  Eirst,  the  acoustic  pressure  is  defined  by 

p{x,  y ,  z)  =  cl)(x,  y  )'F  (z) ,  (20) 
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where  0(x,  j)  is  a  range  dependent  Bessel  funetion,  and  'F(z)  is  the  normal  mode  of  the 
waveguide.  Substitution  into  Equation  3  yields 


0^0  0^0  02'p  ,  , 

'E — -  +  — -  +  0 — -  +  kWo^>  =  0. 


dx^  dy^  0z" 

Separation  of  the  x,  y,  and  z  dependent  terms  yields 


(21) 


1  0^0  1  0^0  1  0"'E  ,2  2  . 

■  + - ^  + - ^  +  k:n  =  0  . 


O  0x"  O  0y2  vp 


(22) 


The  first  and  seeond  term  of  Equation  22  depends  on  x  and  y,  while  the  third  term 
depends  only  on  z.  Therefore,  for  the  equation  to  sum  to  zero  the  first  two  terms  and  the 
last  two  terms  on  the  left  side  must  sum  to  equal  the  same  eonstant  but  of  different  sign. 
Smith  used  the  square  of  the  horizontal  wavenumber  (K)  as  the  separation  variable 
(Smith,  2005),  whieh  yields  the  depth  separated  wave  equation 

0^^^  2  2  7 

^  +  kyw  =  .  (23) 

0Z 


It  is  well  known  that  when  boundaries  are  introdueed  in  layered  media,  the  propagation  is 
eharaeterized  by  the  quantized  modes  defined  by  (Smith,  2005),  whieh  yields 


0z' 


+ky^„,  =Ky„ 


(24) 


Moving  the  right  side  of  the  equation  to  the  left  side  yields  and  faetoring  out  a  kl 
term  yields 


ey 

0z" 


+  kiy-\y, 


y-y. 


=  0 


(25) 


Applying  operators  from  Equation  4  yields 


(26) 
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where  we  ean  introduee  the  operator  q  defined  by  q  =  u  +  s .  Substituting  this  operator 
into  Equation  26  yields 

(27) 

From  Equation  27,  we  ean  deduee  that  a  wave  equation  based  on  the  q  operator  has  the 
normal  mode  shapes  as  its  basis  set  (Smith,  2005).  It  ean  also  easily  be  shown  that  these 
modes  are  identieal  to  the  standard  normal  modes  of  the  Helmholtz  equation  (Smith  and 
Smith,  1998). 

Therefore,  for  sueh  wave  equations  we  ean  assume  the  aeoustie  pressure  field  ean 
be  expressed  as 

N 

p(x,y,z)  =  •  (28) 

n=l 

The  modal  amplitudes  are  then  equal  to  (Smith,  2005) 

"  L  •  (29) 

•"o  p(z) 

The  only  issue  left  is  our  wave  equation,  as  defined  in  Seetion  B,  based  on  the  q 
operator.  From  Equation  13,  the  is  defined  as 

Qop  ~  y/l  +  U  +  V  +  s  .  (30) 

In  a  purely  layered  waveguide,  the  v  operator  ean  be  negleeted  yielding 

Qop  =  'J^  +  U  +  £  ,  (31) 

whieh  would  be  eonsistent  with  the  q  operator.  However,  sinee  the  PE  uses  the 
Thompson  and  Chapman  approximation  the  is  defined  as 

a^«^/^T^+^/^T7-l.  (32) 

Therefore,  the  PE  does  not  use  the  q  operator  exaetly.  Aeeording  to  Smith  and  Smith 
(1998),  this  introduees  seeond  order  errors  but  ean  be  minimized  with  the  ehoiee  of  value 
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for  the  reference  sound  speed  (c^ ).  The  result  is  that  some  modal  amplitudes  may  exhibit 
fluctuations  with  range.  Tests  on  the  stability  of  the  modal  decomposition  were 
performed  in  order  to  minimize  the  impact  of  this  effect  on  the  results  presented  in  this 
thesis. 
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III.  HOW  THE  ANALYTICAL  MODEL  WORKS 


A.  INTRODUCTION 

In  this  chapter,  a  brief  deseription  of  the  Single  Seattering  Analytieal  model  is 
given.  The  deseription  shall  eover  the  pertinent  material  from  an  artiele  Professor  Colosi 
published  in  the  Journal  of  the  Aeoustieal  Soeiety  of  Ameriea  (JASA)  (Colosi,  2008)  and 
from  Chapter  6  of  Introduetion  to  Aeoustie  Fluetuation  (Colosi  and  Flatte,  2008).  The 
first  part  of  the  diseussion  eovers  how  the  expression  for  mode  eoupling  was  developed 
and  how  the  expression  was  altered  to  aeeount  for  multiple  solitons.  Lastly,  an 
explanation  is  given  for  how  the  expression  was  utilized  in  the  model. 


B,  AN  EXPLANATION  OF  HOW  THE  ANALYTICAL  EXPRESSION  WAS 
DEVELOPED  FOR  MODE  COUPLING  INTO  AND  OUT  OF  A  MODE 


In  this  seetion,  an  explanation  is  given  of  how  the  analytieal  expression  for  mode 
eoupling  was  developed.  The  explanation  will  follow  elosely  the  Colosi  JASA  artiele 
previously  refereneed  (Colosi,  2008).  This  referenee  to  the  artiele  is  given  now  to 
prevent  exeessive  referenees  baek  to  the  artiele  in  this  seetion.  The  explanation  is 
reprodueed  again  in  this  paper  instead  of  just  refereneing  the  JASA  artiele  sinee  exeerpts 
will  be  taken  from  other  sourees  in  order  to  provide  a  suffieient  explanation  of  the  theory 
used  in  this  analysis. 


Colosi  used  the  Creamer  formalism  for  two-dimensional  mode  eoupling  in 
shallow-water  to  produee  an  expression  for  the  ehange  in  modal  energy  due  to  aeoustie 
propagation  through  a  soliton  (Creamer,  1996).  The  ssp  was  defined  in  range  (r)  and 
depth  (z)  by 

c{r,z)  =  c{z)  +  dc{r,z) .  (33) 

Here,  c(z)is  the  baekground  ssp  and  5c{r,z)\s  the  perturbation  to  the  ssp  eaused  by  a 
soliton.  The  aeoustie  pressure  field  is  defined  in  terms  of  the  summation  of  the  produet 
of  the  modal  amplitudes  ( (r)  )  and  the  unperturbed  vertieal  mode  shapes  ( (j)^  (z)  ). 


p{r,z) 


fAir) 
«=i  Vv 


(34) 
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Here,  is  the  modal  wave  number.  The  eonsequenee  of  using  the  unperturbed  vertieal 

mode  shape  is  that  all  of  the  aeoustie  variability  is  eonfined  to  the  modal  amplitude 
funetion.  This  faet  will  be  exploited  in  Chapter  IV  of  this  paper  by  allowing  for  a 
eomparison  of  the  analytieal  and  PE  models  by  only  eomparing  the  ehange  in  modal 
energies,  whieh  are  direetly  derived  from  the  modal  amplitudes.  Colosi  uses  the 
quasistatie,  narrow  angle,  weak  forward  seattering  approximation  of  Dozier  and  Tappert 
(1978a),  and  Dozier  (1983)  to  show  that  the  sealed  mode  amplitude  y/^  = 
obeys  the  one-way  equation 

^  (r)  (35) 

dr 

where  /  „  =  /  -  /„  is  the  differenee  between  eomplex  horizontal  wave  numbers.  Note 
that  in  this  eontext,  n  is  used  to  denote  the  mode  of  interest  and  m  is  used  to  denote  the 
mode  to  whieh  n  is  being  eoupled  to.  The  symmetrie  eoupling  matrix  is  defined  as 

Pmnir)  =  kAu{r,z)^^dz.  (36) 

cr(z) 

Here,  cr(z)  is  the  density  profile  as  a  function  of  depth,  u{r,z)  =  Sc{r,z)lc^  is  the 
fractional  sound  speed  variance,  is  the  reference  wave  number,  and  is  the  reference 
sound  speed  (Colosi  and  Flatte,  2008). 

The  change  in  ssp  with  depth  is  typically  small  in  the  natural  environment. 

Therefore,  the  fractional  sound  speed  variance  caused  by  a  soliton,  which  is  a 

displacement  of  the  layer  in  depth,  is  also  typically  small.  Colosi  used  this  fact  to 
produce  a  solution  to  Equation  35  in  terms  of  the  Dyson  Series  (Dyson,  1949;  Sakurai, 
1985).  The  solution  to  second  order  is 

N  N  N 

F.W  =  !<'.(0)-<£F.(0)5l>-££r,(0)5®,  (37) 

m=l  m=l  j=\ 

where  R  corresponds  to  the  receiver  range.  The  first  and  second  order  scattering  matrices 
and  Sll]  are  defined  by 

mn  mnj  J 
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and 


S^mn=\dkp^„{k)\  ^  dr'e 

=  j  dk\  dk'p^„  (k)  *  p.^  (^')  ^'''L/2 


((L.-kyyi-,„-k’)r") 


(38) 

(39) 


For  the  single  scattering  model,  Colosi  concluded  that  the  second  order  scattering  term  in 
Equation  39  is  typically  small  for  the  shallow-water  conditions  under  consideration  and 
therefore  could  be  neglected  without  introducing  a  gross  inaccuracy. 

To  further  develop  the  expression  for  the  mode  coupling  matrix  {p^^{k)),  Colosi 
defined  the  fractional  sound  speed  variance  in  term  of  sound  speed  gradient  {dd  dz  ),  the 
vertical  displacement  of  the  soliton  (i^„),  the  mode  one  vertical  mode  shape  (W^(z)),  the 
parameter  F{r,t)  which  contains  the  variation  in  range  and  time  caused  by  the  soliton. 


uir,z,t)  =  —i^iz)Vw,iz)Fir,t) 


Co  ydz 


Substituting  Equation  40  into  Equation  36  then  yields 

k  dc 


(40) 


Pn 


Xrd)  =  ko\l 


''o  V 


-1  (  Jc  .  .  V  .  . 

UZ  J  CTyZ) 


dz  . 


(41) 


If  the  contribution  to  mode  coupling  from  the  depth  dependence  of  the  soliton  and  the 
mode  structure  is  defined  as 


= 


-^oK 


dz , 


(42) 


c,  ""  j  a(z) 

then  the  mode  coupling  matrix  is  defined  as 

Pmn(Ft)  =  Fir,t)Z^„.  (43) 

Colosi  uses  a  Gaussian  function  to  represent  the  range  and  time  dependence 
associated  with  a  soliton,  where 


F(r,t)  =  exp 


A" 


(44) 


Here,  r^{t)  and  A  are  the  position  and  width  of  the  soliton,  respectfully.  Substitution  of 
Equations  42,  43,  and  44  back  into  Equation  38  and  integrating  yields 
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(45) 


An  expression  for  the  modal  energies  ean  now  be  defined  in  terms  of  the  physieal 
parameters  of  the  soliton.  The  modal  energies  will  be  defined  as  the  square  of  the 
magnitude  of  the  modal  amplitudes.  First,  A„(^)  will  be  substituted  for  in 

Equation  37,  and  then  multiplied  by  its  complex  conjugate.  After  dropping  the  second 
order  terms  for  the  single  scattering  approximation  the  result  is 

N 

=M.(0)|  +2£4(0)4,(0)Ini(5i'’)'  (46) 

m=l 

Substitution  of  Equations  45  into  Equation  46  then  yields 

sin(^^„rj.  (47) 

This  is  the  expression  used  in  the  model  to  determine  the  change  in  modal  energy  after 
acoustic  propagation  through  a  soliton. 

There  are  three  fundamental  assumptions  made  in  the  development  of  this 
theoretical  model  that  are  worth  re-stating  here  as  they  may  impact  the  subsequent 
comparisons  between  the  theoretical  predictions  and  the  solutions  computed  from  the 
3DMMPE  model.  These  assumptions  are: 

1)  the  interactions  of  the  field  with  the  ISW  perturbation  are  only  2-dimensional; 

2)  the  coupling  introduced  by  the  forward  scattering  is  assumed  to  be  relatively 
weak;  and 

3)  only  first  order  coupling  effects  are  considered. 

While  the  second  assumption  may  be  considered  a  reasonable  approximation  of  the 
interactions  of  the  field  with  the  ISW,  the  impact  of  the  first  and  third  assumptions  are 
unknown.  Potential  effects  observed  in  the  analysis  will  be  discussed  later. 


=|4(0)f  +2V^Z4(0)4(0)Z..Af  exp 
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1,  Explanation  of  How  the  Analytical  Expression  Is  Altered  to  Account 
for  Propagation  through  Multiple  Solitons 


In  this  section,  the  terms  are  identified  that  account  for  multiple  solitons.  How 
these  terms  are  altered  and  what  changes  in  the  code  are  required  to  apply  the  alterations. 
The  variables  that  can  change  based  on  the  number  of  solitons  are  ,  A ,  and  r  .  All  of 


these  variables  can  be  accounted  in  our  expression  for  the  modal  energies  by  adding  a 
summation  to  Equation  45  yielding 


S'iH  =  exp 


p=\ 


-k  A 

mn  p 


Jkm„rPp) 


(48) 


Equation  42  changed  by  allowing  the  vertical  displacement  of  the  soliton  to  change  from 
one  soliton  to  another  yielding 


Z 


mn 


-^o(p)K 


cr(z) 


dz . 


(49) 


As  shown  in  Eigure  2,  all  of  the  mode  coupling  effects  from  the  solitons  are  applied  at  the 
axis  of  the  first  soliton  in  the  packet.  Eor  this  reason,  a  line  was  superimposed  on  the 
RME  field  to  identify  where  the  end  of  the  packet  is  located.  As  a  result  of  where  the 
mode  coupling  effects  are  applied,  the  RME  field  is  not  valid  for  the  analytical  model 
until  after  the  soliton  packet.  In  the  analysis  presented  in  this  paper,  comparisons 
between  the  models  are  only  made  for  the  RME  field  after  the  packet. 
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Analytical  Model 
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Figure  2.  Ratio  of  Modal  Energies  (mode  1,  150  Hz,  20“  soliton  angle) 


C.  EXPLANATION  OF  HOW  THE  ANALYTICAL  EXPRESSION  FOR 

MODE  COUPLING  IS  UTILIZED  IN  THE  MODEL 

In  this  section,  the  procedure  used  by  the  analytic  model  to  implement  Equation 
47  is  briefly  explained.  An  explanation  of  how  the  model  deals  with  changing  effective 
soliton  widths  is  also  given.  A  sample  of  the  results  is  also  explained  using  Equation  47. 

Eigure  3  is  an  example  of  a  surface  plot  of  the  RMEs.  Eirst,  the  model  breaks  up 
the  2-D  field  into  angles  of  propagation  from  the  source.  Here,  the  2-D  field  corresponds 
to  the  2-D  surface  plot  that  the  RME  field  will  be  superimposed  onto.  Then,  the  distance 
from  the  source  to  the  soliton  ( )  is  calculated  for  each  angle  of  propagation.  The 
values  are  stored  in  a  vector  that  references  the  distance  to  the  angle  of  propagation.  The 
model  then  inserts  the  physical  parameters  of  the  soliton  and  the  values  into  Equation 

47  and  the  corresponding  value  of  modal  energy  is  stored  in  an  NxM  matrix  referenced  to 
the  angle  of  propagation  and  the  mode  number.  Each  modal  energy  value  is  divided  by 
the  initial  modal  energy  value  at  the  source.  The  resultant  value  is  the  RME  of  the 
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perturbed  and  unperturbed  field.  The  RME  is  then  a  normalized  value,  whieh  gives  the 
amount  of  energy  coupled  into  or  out  of  the  particular  mode. 

Note  that  the  soliton  width  in  Equation  47  corresponds  to  the  effective  soliton 
width  (A^).  Due  to  the  spherical  propagation  the  acoustic  path  will  incident  onto  the 

soliton  at  different  angles.  Therefore,  when  the  soliton  is  oriented  horizontally,  the 
distance  traveled  by  the  acoustic  energy  in  order  to  cross  the  soliton  increases  with 
increasing  angle  of  propagation  from  the  zero  cross-range  axis.  The  model  accounts  for 

the  change  in  the  effective  width  by  ~^cos6’  ^  actual  width  of  the 

soliton  and  6  is  the  angle  measure  between  the  acoustic  path  and  the  vertical  axis. 

The  RME  field  is  divided  up  into  an  NxM  grid  with  each  grid  segment 
corresponding  to  a  cell  within  an  NxM  matrix.  The  distance  from  the  source  to  each  cell 
is  calculated  along  with  angle  of  propagation  from  the  source.  If  r>r^,  then  the  cell  is 

assigned  the  RME  value  corresponding  to  its  angle  of  propagation.  If  r  <r^  then  the  cell 
is  assigned  the  value  of  one.  The  NxM  matrix  is  then  plotted  onto  a  surface  plot  against 
the  down  range  and  cross  range  values  from  the  source  as  in  Eigure  3. 


ANALYTICAL  MOLEL 


Cross  Range  (km) 


Eigure  3.  Ratio  of  Modal  Energies  (mode  1,  150  Hz,  20“  soliton  angle) 
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The  features  in  Figure  3  ean  be  explained  through  the  use  of  Equation  47.  As 
increases  along  the  length  of  the  soliton  with  respect  to  ,  the  second  term  in  the 
expression  oscillates  sinusoidally  with  the  sine  function.  This  explains  the  striation 
pattern  after  the  soliton.  Also,  as  increases  along  the  length  of  the  soliton,  it  causes 

competing  effects  in  the  second  term  of  Equation  47.  Eor  the  range  over  which 

changes  the  dominant  effect  is  to  increase  the  magnitude  of  the  second  term.  Therefore, 
the  magnitude  of  the  ratios  will  increase  with  increasing  A  ^  along  the  length  of  the 

soliton. 
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IV.  VALIDATION  OF  ANALYTICAL  MODEL  WITH  3DMMPE 

MODEL 


A,  COMPARISON  OF  THE  ANALYTICAL  MODEL  WITH  THE  PE 
NUMERICAL  MODEL  FOR  VARYING  ORIENTATIONS  OF  THE 
SOLITON 


In  this  section,  the  results  of  the  PE  numerical  model  and  the  analytical  model 
will  be  compared  for  the  purpose  of  validating  the  predictions  of  the  theory.  The 
comparison  will  be  based  on  the  correlation  of  the  RMEs.  This  approach  is  valid  since 
both  techniques  assume  that  the  acoustic  variability  introduced  by  the  mode  coupling 
effects  is  confined  to  the  modal  amplitudes.  The  RME  is  then  defined  by 


RME  = 


Air) 


perturbed 


Air) 


unperturbed 


(50) 


In  this  notation,  the  unperturbed  modal  energy  refers  to  the  acoustic  field  without  a 
soliton  present  and  perturbed  refers  to  the  acoustic  field  with  a  NEIW  present. 

Two  criteria  will  be  used  to  measure  how  well  the  RMEs  correlate.  The  first  will 
be  a  normalized  correlation  coefficient,  which  is  a  measure  of  how  well  the  different 
regions  of  mode  coupling  align  geometrically  (Rodgers  and  Nicewander,  1988).  The 
second  will  be  a  comparison  of  the  magnitude  of  the  RMEs. 


1,  Comparison  of  the  PE  Numerical  Model  with  the  Analytical  Model 
for  a  Single  Soliton  Perpendicular  to  the  Main  Direction  of  Acoustic 
Propagation 

In  this  section,  the  models  will  be  compared  for  the  case  where  propagation  of  the 
soliton  is  in  the  same  direction  as  the  x-axis.  This  orientation  will  be  referred  to  as  the  0° 
soliton  angle,  where  the  soliton  angle  is  the  angle  between  the  propagation  of  the  soliton 
and  the  central  radial  acoustic  propagation.  In  this  case,  the  central  radial  of  the  acoustic 
energy  will  be  incident  at  a  normal  angle  to  the  soliton,  where  the  mode  coupling  occurs 
entirely  in  the  vertical  plane  (Colosi,  2008).  Eirst,  a  surface  plot  of  the  RME  from  each 
model  will  be  compared  qualitatively  for  geometric  orientation  and  magnitude.  Then,  a 
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line  plot  of  the  RMEs  at  5  km  from  the  souree  will  be  used  in  order  to  make  a 
quantitative  eomparison  of  the  geometrie  orientation  and  magnitude. 

Figure  4  is  the  RME  distribution  for  mode  1  at  150  Hz  produeed  by  the  PE  model. 
Figure  5  is  the  similar  RME  distribution  as  eomputed  by  the  analytieal  model.  Both  plots 
show  a  striation  pattern  in  the  RME  field  after  the  soliton.  The  size  and  general 
orientation  of  the  striation  patterns  are  eomparable.  Speeifieally,  the  striation  patterns  are 
symmetrie  about  the  eenter  of  the  plot  with  roughly  the  same  magnitude  of  RME. 
However,  the  analytieal  model  shows  an  extra  defoeusing  region  starting  at  2  km  eross 
range  just  after  the  soliton.  This  suggests  that  the  width  of  the  striation  patterns  for  the 
PE  is  slightly  wider  than  predieted  by  the  analytieal  model.  The  deviation  from  unity  of 
the  RME  for  the  analytieal  model  is  also  eonsistently  larger  than  for  the  PE.  This  eould 
be  the  result  of  the  model  implementing  the  initial  eonditions  differently.  More  than 
likely,  it  is  a  seeond  order  seattering  effeet  whieh  as  diseussed  in  Chapter  III,  this 
analytieal  model  does  not  aeeount  for.  The  PE  results  in  Figure  4  also  eontain  an 
oseillation  in  the  striation  pattern  that  oseillates  as  a  funetion  of  range  for  a  fixed  angle 
of  propagation.  This  oseillation  is  not  present  in  the  analytieal  predietions.  This  is  likely 
due  to  very  minor  3-D  horizontal  refraetion  effeets,  whieh  are  not  aeeounted  for  in  the  2- 
D  assumption  of  the  theory.  The  eomparison  of  the  surfaee  plots  for  the  other  three 
modes  analyzed  (two,  four,  and  five)  had  similar  qualitative  results. 

It  is  important  to  note  that  Figures  4  and  5  are  the  ratio  of  the  modal  energies 
(RME)  for  the  perturbed  and  unperturbed  aeoustie  fields ||H^(r)|^^^y|H^(r)| 

Therefore,  this  is  a  eomparison  of  the  amount  of  mode  eoupling  that  oeeurs  in  eaeh 
model,  not  a  direet  eomparison  of  the  modal  amplitude  distribution  of  eaeh  model.  Also, 
the  striation  pattern  present  in  eaeh  figure  is  made  up  of  foeusing  and  defoeusing  regions 
as  a  result  of  vertieal  mode  eoupling  eaused  by  propagation  through  the  soliton.  The 
physieal  eonsequenee  of  the  sinusoidal  fluetuation  is  that  the  RME  field  is  made  up  of 
regions  where  modal  energy  is  transferred  to  or  from  one  mode  to  another.  Modal  energy 
is  transferred  to  the  mode  in  the  foeusing  regions  and  out  of  the  mode  in  the  defoeusing 
region. 
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PE  Model 
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Figure  4.  Ratio  of  Modal  Energies  (mode  1,  150  Hz,  0°  soliton  angle) 
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Figure  5.  Ratio  of  Modal  Energies  (mode  1,  150  Hz,  0°  soliton  angle) 
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Some  explanation  should  be  given  for  the  oscillation  in  the  focusing  pattern 
shown  in  Figure  4,  and  why  it  does  not  show  up  in  Figure  5.  For  the  Single  Scattering 
Approximation  used  in  the  analytical  model,  all  of  the  mode  coupling  effects  of  the 
soliton  is  accounted  for  at  the  center  of  the  soliton  located  at  2  km  down  range  in  Figure 
4.  Therefore,  the  analytical  model  would  not  account  for  coupling  effects  down  range  of 
the  soliton  without  accounting  for  second  order  or  higher  scattering  effects  (Colosi, 
2008).  The  fluctuations  in  Figure  4  for  the  PE  model  are  possibly  a  numerical  artifact, 
since  they  occur  in  regions  where  the  ssp  is  not  a  function  of  range  from  the  source. 
However,  such  fluctuations  are  not  present  at  higher  soliton  angles  (see,  for  example. 
Figure  11)  and  suggest  that  this  feature  has  a  physical  manifestation.  It  could  be  a  3-D, 
horizontal  refraction  effect  which  the  2-D  model  would  not  compute.  The  PE  model  was 
not  run  in  an  Nx2-D  implementation  to  confirm  this,  unfortunately.  Euture  work  will 
attempt  to  address  this  effect,  both  numerically  and  theoretically.  The  effect  is  not 
significant  for  modes  1  and  2,  however,  in  modes  4  and  5  the  oscillations  are  closer 
together  and  change  the  look  of  the  field  considerably.  Eigures  6  and  7  are  the  surface 
plots  for  mode  5  produced  by  the  PE  model  and  analytical  model. 


PE  Model 


-  1.2 


0.8 


Eigure  6. 


Ratio  of  Modal  Energies  (mode  5,  150  Hz,  0°  soliton  angle). 
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Analytical  Model 
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Figure  7.  Ratio  of  Modal  Energies  (mode  5,  150  Hz,0“  soliton  angle). 


A  quantitative  comparison  of  the  correlations  can  be  made  by  looking  only  at  the 
RME  for  both  predictions  at  a  specific  range.  Eigure  8  is  a  line  plot  of  the  RMEs  at  5  km 
down  range,  with  the  PE  ratio  in  blue  and  the  analytical  ratio  in  red.  The  range  for  the 
comparison  was  chosen  to  be  at  5  km  to  stay  consistent  with  other  line  comparisons 
presented  later  in  this  paper.  A  range  of  3  km  may  have  been  a  better  choice  for 
comparing  the  results  for  this  soliton  angle.  However,  the  results  did  not  change 
appreciably  when  analyzed  at  3  km. 

As  noted  from  the  surface  plots,  Eigure  8  shows  that  the  predicted  amount  of 
mode  coupling  is  consistently  greater  for  the  analytical  model.  The  largest  difference  in 
the  magnitude  of  mode  coupling  occurred  at  0  km  cross  range  and  resulted  in  8%  more 
mode  coupling  in  the  analytical  model.  The  regions  of  mode  coupling  geometrically 
align  very  well  with  a  correlation  coefficient  of  0.9997. 
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Figure  8.  Line  Comparison  of  Ratio  of  Modal  Energies  at  5km  (mode  1,  150  Hz,  0° 

soliton  angle). 


For  every  mode  analyzed  at  this  soliton  angle,  the  magnitude  of  ratios  deviation 
from  unity  is  larger  for  the  analytical  model  than  for  the  PE  model.  This  indicates  that 
the  analytical  model  consistently  predicts  a  larger  amount  of  mode  coupling.  Again,  it  is 
unclear  if  this  is  due  to  limiting  assumptions  of  the  theory  or  numerical  issues  with  the  PE 
model.  The  percent  that  the  RME  deviation  from  unity  is  larger  for  the  analytical  model 
than  for  the  PE  model  for  each  mode  at  a  0°  soliton  angle  are  tabulated  in  Table  1.  Of 
note  is  the  value  for  mode  5,  which  is  almost  twice  as  large  as  its  nearest  value  8%  for 
mode  1.  Table  2  contains  the  correlation  coefficients  at  each  mode  analyzed  at  a  0° 
soliton  angle.  The  data  compiled  in  Tables  1  and  2  indicates  that  the  predictions 
produced  by  the  analytical  model  are  reasonably  consistent  with  the  PE  model  for  the  0° 
soliton  angle  case. 
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Soliton 

Angle 

Mode  1 

Percent  Larger 

Mode  2  Mode  4 

Mode  5 

0 

7.99 

6.57  3.61 

15.27 

Table  1 .  Maximum  Percent  that  the  Analytical  RME  Deviation  from  Unity  is 

Greater  than  the  PE  RME  Deviation  Erom  Unity 


Soliton 

Angle 

Correlation  Coefficient 

Mode  1  Mode  2  Mode  4  Mode  5 

0 

0.9997 

0.9757 

0.9347 

0.8165 

Table  2.  Correlation  Coefficients  for  Perturbed/Unperturbed  Einear  Ratios  at  5  km. 

2,  Comparison  of  Analytical  Model  with  the  PE  Model  for  a  Soliton  at 
Varying  Orientations  to  the  Main  Direction  of  Acoustic  Propagation 

A  surface  plot  of  the  RME  is  used  in  this  section  to  compare  the  analytical  model 
with  the  PE  model  predictions  at  varying  soliton  angles.  Eigures  9  and  10  are  the  RMEs 
for  mode  1,  150  Hz,  at  a  20“  soliton  angle  produced  from  the  PE  model  and  the  analytical 
model,  respectively.  Just  as  in  the  case  for  the  0“  soliton  angle,  the  20“  soliton  angle 
causes  a  striation  pattern  made  up  of  focusing  and  defocusing  regions  in  the  RME  field 
after  the  soliton.  The  patterns  are  due  to  vertical  mode  coupling  caused  by  the  soliton.  In 
fact,  the  work  of  Roush  (2008)  showed  that  horizon  refraction  effects  do  not  become 
significant  until  soliton  angles  of  approximately  85“.  Therefore,  since  the  largest  soliton 
angle  analyzed  in  this  paper  is  80“,  all  of  the  mode  coupling  effects  in  this  paper  will  be 
the  result  of  vertical  mode  coupling. 

The  surface  plots  in  Eigures  9  and  10  correlate  well  over  the  negative  values  of 
cross  range.  The  extra  focusing  region  from  the  analytical  predictions,  displayed  in 
Eigure  10,  indicates  that  the  striation  patterns  are  slightly  wider  as  computed  by  the  PE 
model,  displayed  in  Eigure  9.  The  oscillation  in  Eigure  9  does  not  appreciably  change  the 
RME  field,  and  therefore  is  not  a  major  difference  between  the  models  for  this  soliton 
angle  and  mode.  The  correlation  breaks  down  for  the  positive  cross  range  values.  In 
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Figure  9,  the  PE  model  predicts  that  a  defocusing  region  covers  the  entire  area,  where  the 
analytical  model  in  Figure  10  exhibits  an  extra  focusing  region. 


PE  Model 


Cross  Range  (km) 


Figure  9.  Ratio  of  Modal  Energies  (mode  1,  150  Hz,  20°  soliton  angle) 

ANALYTICAL  MOLEL 


Cross  Range  (km) 


Eigure  10.  Ratio  of  Modal  Energies  (mode  1,  150  Hz,  20°  soliton  angle) 
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Next,  the  50°  soliton  angle  ease  is  examined.  Figures  11  and  12  are  the  RME 
distributions  from  the  PE  model  and  the  analytieal  model,  respeetively,  for  mode  1,  150 
Hz,  and  a  50°  soliton  angle.  A  comparison  between  Eigures  1 1  and  12  does  show  that  the 
models  continue  to  correlate  well.  The  size  of  the  striation  patterns  actually  correlate 
better  over  the  entire  spectrum  of  cross  range  values.  However,  the  difference  in  the 
magnitude  of  the  RMEs  deviation  from  unity  is  more  pronounced  for  this  soliton  angle. 
The  analytical  model  still  shows  a  defocusing  region  in  the  positive  cross  range  area  that 
is  not  present  for  the  PE  model,  but  this  region  is  beyond  the  60°  boundary  condition  of 
the  PE.  The  PE  does  not  produce  accurate  results  at  an  angle  greater  than  60°  from  the 
vertical  due  to  limited  angular  extent  of  the  source.  The  oscillation  in  the  focusing  region 
of  Eigure  1 1  is  more  pronounced  for  this  angle  than  for  the  20°  soliton  angle  but  is  not  as 
significant  to  the  overall  look  of  the  field  as  for  the  0°  soliton  angle.  The  oscillation  is 
also  occurring  in  different  areas  of  the  RME  field  for  each  soliton  angle.  All  of  the  above 
considered,  the  models  still  correlate  well. 

Overall,  from  looking  at  the  surface  plots  of  the  RMEs  presented  here,  the 
analytical  and  PE  numerical  predictions  appear  to  correlate  well.  The  two  models  are  not 
expected  to  produce  point  wise  accurate  results.  Instead,  the  goal  is  to  confirm  that  the 
two  methods  produce  the  same  general  shape  and  magnitude  of  RME.  From  the  above 
surface  plots,  this  goal  appears  to  have  been  reached.  Surface  plots  of  the  RME  at  10° 
soliton  angle  increments  from  0°  to  80°  at  mode  1,  2,  4,  and  5  were  evaluated  as  part  of 
this  analysis.  There,  results  were  consistent  with  the  conclusions  presented  for  the 
surface  plots  shown  above.  One  difference  of  note  is  that  the  oscillation  contained  within 
the  striation  patterns  for  the  PE  model  became  more  significant  in  the  appearance  of  the 
RME  field  with  increasing  mode  number. 
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Figure  11.  Ratio  of  Modal  Energies  (mode  1,  150  Hz,  50“  soliton  angle) 


Analytical  Model 
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Figure  12.  Ratio  of  Modal  Energies  (mode  1,  150  Hz,  50“  soliton  angle) 
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As  for  the  zero  soliton  angle  ease,  a  quantitative  look  at  the  RME  is  taken  by 
looking  at  their  values  for  a  fixed  down  range  value.  Figures  13  and  14  show  how  well 
the  PE  model  and  analytieal  model  eorrelate  for  the  20°  and  50°  soliton  angles.  The 
correlation  ratios  are  0.9893  and  0.8758,  respectively.  However,  a  trend  can  be  seen  in 
the  line  comparisons.  As  the  soliton  angle  increases,  the  difference  in  the  magnitude  of 
the  RMEs  deviation  from  unity  increases.  The  difference  in  the  magnitudes  is  most 
likely  due  to  two  possible  causes.  There  may  be  a  process  that  the  analytical  model  does 
not  account  for  that  is  causing  less  vertical  mode  coupling  to  occur  than  expected.  Or, 
the  method  in  which  the  analytical  model  and  the  PE  model  account  for  the  initial 
conditions  is  different. 


Figure  13.  Fine  Comparison  of  the  Ratio  of  Modal  Energies  at  5km  (mode  1,  150  Hz, 

20°  soliton  angle). 
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Figure  14.  Line  Comparison  of  the  Ratio  of  Modal  Energies  at  5km  (mode  1,  150  Hz, 

50°  soliton  angle). 

Table  3  lists  all  of  the  correlation  ratios  for  the  varying  soliton  angles  for  modes  1, 
2,  4,  and  5.  Mode  3  was  not  used  for  comparison  due  to  the  small  amount  of  acoustic 
energy  initially  put  into  that  mode.  With  the  exception  of  mode  4  at  a  soliton  angle  of  80°, 
all  of  the  correlation  values  are  greater  than  75%.  This  value  can  be  explained  by  referring  to 
Figure  15.  The  large  blue  spike  is  the  PE  RME  between  -1  and  -0.5  km  cross  range  and  is 
caused  by  the  presence  of  the  soliton.  When  the  PE  calculates  the  RME,  it  uses  the 
unperturbed  vertical  mode  shapes.  However,  along  the  soliton  crest,  the  mode  shapes  are 
distorted  due  to  the  change  in  sound  speed  profde.  The  PE  model  assumes  the  same  mode 
shape  at  the  soliton  to  determine  the  RME  due  to  the  complex  calculation  that  would  be 
required  to  calculate  the  exact  mode  shape  at  the  soliton.  This  assumption  is  justified 
considering  the  small  width  of  the  soliton  (100m)  compared  to  the  range  scale  being  used 
in  the  model  (5  km)  and  the  reduced  run  time  required  to  run  the  PE.  However,  this 
assumption  is  the  cause  for  the  divergence  between  the  PE  model  and  analytical  model 
right  at  the  soliton  (i.e.,  the  large  blue  spike).  This  means  that  the  models  do  agree  to  a 
high  degree  of  certainty  in  the  geometric  location  of  the  vertical  mode  coupling. 
However,  the  modal  decomposition  of  the  PE  field  is  not  strictly  valid  in  the  vicinity  of 
the  soliton. 
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Soliton 

Correlation  Coefficient 

Angle 

Mode  1 

Mode  2 

Mode  4 

Mode  5 

0 

0.9997 

0.9757 

0.9347 

0.8165 

10 

0.9810 

0.9810 

0.9046 

0.9128 

20 

0.9894 

0.9831 

0.8183 

0.8850 

30 

0.9715 

0.9652 

0.8712 

0.8294 

40 

0.9509 

0.9647 

0.8684 

0.8107 

50 

0.8758 

0.9281 

0.8155 

0.7499 

60 

0.8191 

0.9103 

0.8187 

0.8417 

70 

0.8552 

0.9331 

0.7985 

0.8262 

80 

0.8438 

0.8028 

0.2828 

0.7671 

Table  3.  Correlation  Coefficients  for  Perturbed/Unperturbed  Linear  Ratios  at  5  km 


Figure  15.  Line  Comparison  of  the  Ratio  of  Modal  Energies  at  5km  (mode  4,  150  Hz, 

80°  soliton  angle). 

For  the  20°  soliton  angle,  the  magnitude  of  the  deviation  from  unity  for  the 
analytical  model  is  7.52%  larger  than  for  the  PE  model.  As  can  be  seen  in  Figures  8,  13, 
and  14  for  mode  1,  there  is  a  good  correlation  between  the  magnitudes  of  the  ratios. 
Table  4  shows  how  much  larger  the  magnitude  of  the  RME  is  for  the  analytical  model 
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than  for  the  PE  model  at  the  maximum  difference.  Overall,  Table  4  shows  an  increase  in 
the  difference  between  the  magnitudes  for  increasing  soliton  angles.  Mode  5  also  has 
larger  differences  in  magnitude  than  the  other  modes.  This  could  be  a  result  of  the  low 
initial  energy  deposited  into  that  mode  due  to  the  source  depth.  Figures  13  and  14  are 
good  examples  of  the  average  mismatch  in  magnitude  of  the  RMEs  with  values  of  7.52% 
and  15.92%.  Figure  16  is  an  example  of  the  more  extreme  mismatch,  with  the  analytical 
model  RME  deviation  from  unity  more  than  40%  larger  than  the  PE  model  RME 
deviation  from  unity. 


Soliton 

Percent  Larger 

Angle 

Mode  1 

Mode  2 

Mode  4 

Mode  5 

0 

7.99 

6.57 

3.61 

15.27 

10 

8.31 

9.88 

6.60 

9.93 

20 

7.52 

13.77 

8.99 

26.70 

30 

7.58 

15.02 

8.19 

24.08 

40 

9.12 

11.43 

8.28 

41.82 

50 

15.92 

11.79 

6.62 

41.99 

60 

18.73 

11.67 

7.81 

34.19 

70 

14.51 

24.99 

7.73 

17.92 

80 

12.52 

14.23 

2.34 

24.95 

Table  4.  Maximum  Percent  that  the  Analytical  RME  Deviation  from  Einity  is 

Greater  than  the  PE  RME  Deviation  from  Einity 
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Figure  16.  Line  Comparison  of  the  Ratio  of  Modal  Energies  at  5km  (mode  5,  150  Hz, 

50  soliton  angle). 

3,  Conclusions 

In  this  seetion,  we  have  analyzed  the  correlation  between  the  RMEs  produced 
from  the  results  of  the  analytical  and  PE  models  for  sound  propagation  through  a  single 
soliton.  Through  those  analyses,  we  have  shown  that  the  analytical  model  does  produce 
results  consistent  with  the  PE  model  over  the  soliton  angles  and  mode  numbers  analyzed. 
Eor  example,  both  models  produce  RME  fields  with  striation  pattern  of  comparable  size, 
shape,  and  magnitude.  However,  some  limitations  have  also  been  discovered.  Eor 
certain  mode  numbers  and  soliton  angles,  the  RME  deviations  from  unity  are  more  than 
49%  larger  for  the  analytical  model  than  for  PE  model.  The  correlation  coefficients  for 
the  geometric  alignment  of  the  striation  pattern  also  dropped  to  as  low  as  0.7671.  These 
differences  in  the  results  do  not  have  a  discernible  trend  in  which  to  predict  conditions 
where  large  errors  would  occur.  The  PE  model  also  contains  patterns  in  the  striation 
patterns  that  oscillate  as  a  function  of  range  from  the  source  for  a  given  angle  of  acoustic 
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propagation  that  did  not  occur  in  the  analytical  model.  These  limitations  make  point  wise 
agreement  of  solutions  unlikely.  However,  the  analytical  model  is  sufficient  at  predicting 
the  same  trends  and  general  behavior  as  produced  by  the  PE  model. 

B,  COMPARISON  OF  ANALYTICAL  MODEL  WITH  THE  PE  MODEL 
PREDICTION  FOR  MODE  COUPLING  RESONANCE  BEHAVIOR 

In  this  section,  to  determine  if  the  PE  model  follows  the  mode  coupling  resonance 
behavior  called  for  by  the  theory,  a  comparison  is  made  between  the  PE  model  and  the 
analytical  model.  Eirst,  the  theory  behind  why  the  resonance  occurs  will  be  reviewed. 
Then,  the  conditions  for  the  comparison  will  be  explained.  Einally,  the  results  from  the 
PE  model  will  be  analyzed  to  determine  how  well  they  compare  with  the  analytieal 
prediction. 


1,  Theory  Behind  the  Mode  Coupling  Resonance  Behavior 


According  to  Colosi  (2008)  the  vertical  mode  coupling  should  exhibit  a  resonance 
eondition  dependent  on  the  soliton  width  (A)  and  the  frequeney  of  the  acoustic  wave  (m). 
Colosi  (2008)  showed  that  the  modal  energy  after  propagation  through  a  soliton  is 

^mn^f 

4 

I  |2 

where  |^„(0)|  is  the  initial  mode  energy,  Amn  is  the  initial  mode  amplitude  matrix,  Z^n  is 

the  matrix  that  incorporates  the  depth  structure  of  the  soliton  for  modes  m  and  n,  kmn  is 
the  beat  wavenumber,  and  ro  is  the  distance  from  the  source  to  the  soliton.  At  small  and 
large  values  of  A,  the  second  term  on  the  right  hand  side  of  Equation  51  is  driven  to  zero. 

Therefore,  a  resonance  condition  is  reached  when  A  -  yfl  I  Work  done  by  Colosi 


and  Elatte  (2008)  showed  that  the  RMS  variation  in  mode  energy  caused  by  one  soliton, 
averaged  over  a  span  of  ro  values,  per  meter  of  internal  displacement,  is  represented  by 


r„(®,A)  =  - 


fAMz  A  exp  - 


^mn^f 


Note  that  here  A^is  used  to  denote  the  effective  width  of  the  soliton  in  the  direction  of 


acoustic  propagation  other  than  normal  to  the  soliton.  This  expression  was  used  to 
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produce  Figure  17,  whieh  was  used  to  determine  the  approximate  width  of  the  soliton 
when  the  resonanee  oeeurs.  Note  that  Figure  17  is  averaged  over  a  set  of  values.  Also, 

the  water  depth  used  to  produee  Figure  14  was  100  m  whieh  is  the  same  water  depth  used 
by  the  PE,  but  the  souree  is  loeated  at  90  m  rather  than  the  50  m  used  by  the  PE. 
Therefore,  the  A  determined  from  Eigure  17  is  not  expeeted  to  yield  exaet  results  for  A 
at  the  resonanee  but  will  provide  a  good  estimate  that  ean  be  used  to  determine  if  the  PE 
is  eonsistent  with  the  analytieal  model. 


Mode  1 


Mode  2 


Frequency  (Hz) 


Eigure  17.  Contours  of  RMS  Relative  Mode  Energy  Eluetuation  per  Meter  of  Soliton 
Displaeement,  rn(co,A)  (m'^)  for  a  100  m  Water  Depth  Canonieal  Shallow  Water 
Oeean  Waveguide  as  Deseribed  by  Colosi  and  Elatte  (2008). 


35 


2,  PE  Mode  Coupling  Resonance 

The  theory  suggests  that  Figure  17  above  can  be  used  to  determine  the 
approximate  location  of  the  vertical  mode  coupling  resonance  based  on  the  frequency  and 
width  of  the  soliton.  Therefore,  the  PE  model  can  be  used  to  validate  the  analytical 
model  by  varying  the  soliton  width  and/or  acoustic  frequency  to  determine  if  the 
resonance  occurs,  as  suggested  by  the  theory.  In  this  analysis,  a  constant  acoustic  source 
frequency  of  150  Hz  is  used,  and  the  soliton  width  in  the  PE  model  (A^^)  is  varied 

between  50  m  and  500  m  in  increments  of  50  m.  Note  that  the  results  for  the  200  m 

are  not  used  due  to  a  processing  problem  that  could  not  be  resolved.  Erom  Eigure  17 
above,  the  vertical  mode  coupling  resonance  for  an  acoustic  signal  at  150  Hz  should 
occur  at  an  approximate  soliton  width  of  300m  for  mode  1  and  250m  for  mode  2.  Eigure 
18  shows  the  RMEs  at  a  fixed  down  range  value  of  5  km  for  varying  widths.  All  of  the 
solitons  are  oriented  at  a  soliton  angle  of  45°.  A  similar  plot  is  shown  in  Eigure  19  for 
mode  2.  Erom  Eigure  18,  it  is  apparent  that  a  resonance  condition  does  exist.  In  fact, 
one  resonance  condition  exists  for  each  focusing  and  defocusing  region  in  both  modes. 
The  A^^  at  resonance  for  the  defocusing  region  at  1  km  cross  range  in  Eigure  18  is  250 
m.  That  is  less  than  a  17%  error  from  the  300  m  width  predicted  from  Eigure  17. 

Note  from  Eigures  18  and  19  that  the  values  of  the  Ap^  at  the  resonances  are  not 
the  same  for  each  focusing  and  defocusing  region.  The  different  Ap^  at  resonance  is 

caused  by  the  spherical  wave  fronts  intersecting  the  soliton  at  different  angles  along  its 
length.  Eor  this  analysis,  the  acoustic  source  is  a  point  source  centered  at  0  km  cross 
range  and  0  km  down  range.  The  soliton  is  oriented  at  a  45°  angle  from  the  horizontal 
and  centered  at  a  point  located  at  0  km  cross  range  and  2  km  down  range.  The  angle 
between  the  path  the  acoustic  energy  takes  and  the  soliton  changes  along  the  length  of  the 
soliton.  Therefore,  the  distance  the  acoustic  energy  travels  while  crossing  through  the 
soliton  changes  along  the  length  of  the  soliton. 

Due  to  the  symmetry  along  its  axis,  the  change  in  distance  through  the  soliton  can 
be  accounted  for  with  an  effective  soliton  width  (A^  =  A/cos(6*^))  (Colosi,  2007).  The 
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resonance  at  1  km  cross  range  is  oriented  such  that  the  angle  between  the  acoustic  path 
and  the  normal  to  the  soliton  is  zero.  Therefore,  the  is  equal  to  the  .  For  the 

remaining  resonances,  since  the  soliton  width  read  from  Figure  17  is  the  A^  and  the 

width  read  from  Figures  18  and  19  are  the  actual  widths  of  the  soliton,  to  make  a 
comparison  the  actual  widths  of  the  soliton  must  be  converted  to  effective  widths.  For 
the  focusing  region  located  at  about  -750  m  cross  range,  the  angle  between  the  path  of 
propagation  and  the  normal  to  the  soliton  is  about  60°.  Therefore,  the  effective  soliton 
width  is 

A^  =  A/cos(6’^)  =  150m/cos(60.26)  =  302m.  (53) 

The  percent  difference  from  the  300  m  that  was  predicted  from  Figure  17  is  less 
than  one  percent.  The  effective  soliton  widths  for  all  of  the  resonance  shown  in  Figures 
18  and  19  are  tabulated  in  Table  5.  All  of  the  effective  soliton  widths  are  within  20%  of 
the  predicted  soliton  width  except  two,  which  are  located  at  -1.7  km  and  -2.25  km  cross 
range  in  mode  2.  The  percent  errors  for  these  soliton  widths  are  about  70  and  40  percent, 
respectfully.  Overall,  the  results  do  show  that  the  PE  model  produces  results  that  are 
consistent  with  the  analytical  model. 
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Figure  18.  Comparison  of  Mode  Coupling  for  Different  Soliton  Widths  at  5  km 

(mode  1,  150  Hz,  45  soliton  angle) 
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Figure  19.  Comparison  of  Mode  Coupling  for  Different  Soliton  Widths  at  5  km(mode 

2,  150  Hz,  45  soliton  angle) 


Mode  1 

Mode  2 

Cross  Range 

Ape 

Aeff 

%  Diff 

Cross  Range 

Ape 

Aeff 

%  Diff 

1  km 

250 

250 

16.67% 

1 .4  km 

275 

296 

18.40% 

-0.75  km 

150 

302 

0.67% 

-0.75  km 

150 

296 

18.40% 

-1.65  km 

125 

348 

16.00% 

-1.7  km 

150 

424 

69.60% 

-2.2  km 

100 

335 

11.67% 

-2.25  km 

100 

335 

34.00% 

Table  5.  Effective  Soliton  Width  for  each  Resonance  Shown  in  Figure  18  and  19 


3,  Conclusions 

In  this  section,  the  mode  coupling  resonance  condition  for  the  PE  was  compared 
to  the  analytical  predictions  for  varying  soliton  widths.  The  results  tabulated  in  Table  5 
show  that  the  PE  agrees  well  with  the  theory  for  mode  1.  However,  the  PE  begins  to 
differ  considerably  from  the  analytically  predicted  values  for  mode  2.  The  fact  that  one 
mode  agrees  and  the  other  mode  does  not  is  not  enough  data  to  really  determine  how  well 
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the  models  correlate  as  a  whole.  Further  analysis  should  be  conducted  to  include  more 
modes  and  reproduce  Figure  17  using  the  correct  initial  conditions. 


C.  COMPARISON  OF  THE  ANALYTICAL  MODEL  WITH  THE  PE  MODEL 

FOR  MULTIPLE  SOLITONS  AT  VARYING  ORIENTATIONS 

In  this  section,  the  results  of  the  PE  model  and  analytical  model  will  be  compared 
for  the  purpose  of  validating  the  results  of  the  analytical  model  for  a  soliton  packet.  As  in 
Section  A,  the  comparison  will  be  based  on  how  well  the  RME  from  each  model  correlate 
with  the  other.  This  approach  is  still  valid  for  soliton  packets,  since  both  models  still 
assume  that  the  acoustic  variability  introduced  by  the  soliton  packet  is  confined  to  the 
RME. 

The  same  criteria  will  be  used  to  measure  the  correlation  between  models  for  the 
soliton  packet  case  as  for  the  single  soliton  case.  The  first  criteria  will  be  a  correlation 
coefficient,  which  is  a  measure  of  how  well  the  different  regions  of  mode  coupling  align 
geometrically  (Rodgers  and  Nicewander,  1988).  The  second  criteria  will  be  a 
comparison  of  the  magnitude  of  the  RME  at  the  maximum  difference  between  them. 

1.  Comparison  of  the  Analytical  Model  with  the  PE  Model  for  a  Soliton 
Packet  Perpendicular  to  the  Main  Direction  of  Acoustic  Propagation 

In  this  section,  the  models  will  be  compared  for  the  case  where  the  soliton  packet 
is  oriented  at  a  0“  soliton  angle.  As  noted  in  Section  A,  the  soliton  angle  is  measured 
from  the  horizontal  plane  which  is  perpendicular  to  the  main  direction  of  acoustic 
propagation.  The  soliton  packet  is  composed  of  three  solitons,  each  100  m  in  width  and 
separated  by  500  m.  The  vertical  displacement  of  the  solitons  from  the  first  soliton  to  the 
third  is  10  m,  9  m,  and  8  m,  respectfully.  In  this  notation,  a  positive  vertical  displacement 
of  the  solitons  is  in  the  positive  z-direction,  which  is  positive  with  increasing  depth 
(down  direction).  The  0“  soliton  angle  case  is  ideally  the  simplest  case  since  the  majority 
of  the  acoustic  energy  will  be  incident  at  a  near  normal  angle  to  the  soliton  where  the 
mode  coupling  occurs  entirely  in  the  vertical  plane  (Colosi,  2008).  As  in  the  single 
soliton  case  a  surface  plot  of  the  RME  for  each  model  will  be  compared  qualitatively  for 
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geometric  orientation  and  magnitude.  A  line  plot  of  the  RME  at  5  km  from  the  source 
will  be  used  to  make  a  quantitative  comparison  of  the  geometric  orientation  and 
magnitude  of  the  RME. 

Eigures  20  and  21  are  surface  plots  of  the  RMEs  produced  by  the  PE  model  and 
the  analytical  model.  In  Figure  20,  the  end  of  the  soliton  packet  is  easily  identified  by  the 
change  in  the  magnitude  of  the  RME  at  2  km  down  range  from  the  source.  For  the 
analytical  predictions,  all  of  the  effects  of  the  soliton  packet  are  summed  and  applied  at 
the  leading  edge  of  the  soliton  packet  (Colosi,  2008).  Therefore,  the  RME  field  produced 
by  the  analytical  model  is  not  valid  until  after  the  final  soliton  in  the  packet.  To  identify 
this  point,  a  line  has  been  superimposed  onto  the  RME  field  in  Figure  21  to  identify  the 
end  of  the  soliton  packet.  The  PE  model  also  has  areas  where  the  RME  field  is  not  valid. 
At  the  interface  of  each  individual  soliton  the  PE  assumes  that  the  same  unperturbed 
vertical  mode  shapes  exist  across  the  solitons.  This  is  not  the  case,  since  the  soliton 
distorts  the  mode  shapes  across  their  interface.  Therefore,  the  RME  are  not  valid  in  parts 
of  the  field  where  a  soliton  is  located.  However,  this  assumption  is  valid  over  most  of  the 
field  ,due  to  the  small  area  an  individual  soliton  covers,  compared  to  the  rest  of  the  field 
and  the  faster  runtime  the  program  gains  by  not  attempting  to  recalculate  the  mode 
shapes.  For  the  purpose  of  this  analysis,  the  models  will  only  be  compared  for  the  part  of 
the  field  after  the  last  soliton  in  the  packet. 

The  RME  fields  produced  by  the  PE  model  in  Figure  20  and  the  analytical  model 
in  Figure  21  have  the  same  general  shape.  Both  fields  consist  of  striation  pattern 
symmetric  about  the  center.  The  oscillation  in  the  focusing  region  of  Figure  20  does 
significantly  affect  the  appearance  of  the  field.  The  width  of  the  focusing  regions  in 
Figure  20  also  does  not  appear  to  be  as  wide  as  the  corresponding  regions  in  Figure  21. 
This  is  apparent  by  the  fact  that  Figure  20  has  two  focusing  regions  per  side  but  Figure  21 
only  has  one.  Also,  the  defocusing  regions  in  Figure  21  at  about  ±2.5  km  cross  range  is 
much  more  intense  than  in  Figure  20. 


40 


PE  Model 


Cross  Range  (km) 


Figure  20.  Ratio  of  Modal  Energies  (mode  1,  150  Flz,  0°  soliton  angle) 


Analytical  Model 


Down  Range  (km) 


Figure  21.  Ratio  of  Modal  Energies  (mode  1,  150  Hz,  0°  soliton  angle) 
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A  comparison  of  the  PE  model  and  analytical  model  for  surface  plots 
eorresponding  to  modes  2,  4,  and  5  were  performed  with  similar  results.  However,  the 
analysis  does  yield  two  major  differences  between  the  models.  The  first  difference  is  that 
the  RME  from  the  PE  model  ean  vary  greatly  from  the  analytical  model  as  a  result  of  an 
oscillation  within  the  striation  pattern.  The  analytical  model  does  not  contain  these 
oseillations.  Eigures  22  and  23  demonstrate  how  signifieantly  an  effeet  these  oseillations 
can  have  on  the  RME  field.  The  seeond  differenee  is  that  the  analytical  model  usually 
has  a  larger  amount  of  mode  eoupling  than  the  PE  model.  This  is  demonstrated  by  the 
larger  magnitude  of  deviation  form  unity  in  the  RME  field  for  the  analytieal  model.  Both 
differences  will  be  addressed  again  later  in  this  paper  using  the  correlation  eoefficient  and 
by  comparing  the  maximum  differenee  in  model  amplitudes. 


PE  Model 


Cross  Range  (km) 


Eigure  22.  Ratio  of  Modal  Energies  (mode  5,  150  Hz,  0“  soliton  angle) 
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Figure  23.  Ratio  of  Modal  Energies  (mode  5,  150  Flz,  0°  soliton  angle) 


A  quantitative  comparison  of  the  correlations  shall  now  be  made  by  looking  only 
at  the  RMEs  at  a  specific  down  range  value.  Eigure  24  is  a  line  plot  of  the  RMEs  with 
the  PE  model  in  blue  and  the  analytical  predictions  in  red.  The  range  for  the  comparison 
was  chosen  to  be  at  5  km  to  stay  consistent  with  other  line  comparisons  presented 
elsewhere  in  the  paper.  The  RMEs  in  Eigure  24  align  very  closely,  with  a  small 
divergence  at  cross  ranges  greater  than  2  km.  At  the  maximum  difference  between 
RMEs,  the  RME  deviation  from  unity  is  only  1.55%  greater  for  the  analytical  model  than 
for  the  PE  model.  The  RMEs  for  this  mode  have  a  correlation  coefficient  of  0.9881. 
Therefore,  for  this  mode  computed  with  this  soliton  angle,  at  this  rang,  the  analytical 
model  produces  results  very  consistent  with  the  PE  model.  Tables  6  and  7  contain  the 
tabulated  maximum  percent  difference  and  correlation  coefficients  for  modes  1,  2,  4  and 
5  at  the  0”  soliton  angle.  Table  7  shows  that  the  geometric  alignment  of  the  mode 
coupling  regions  decreases  with  increasing  mode  number.  Similarly,  Table  6  shows  that 
the  magnitude  of  the  RMEs  deviation  from  unity  increases  with  increasing  mode  number. 
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Figure  25  is  the  line  eomparison  for  mode  5  and  is  the  case  for  the  worst  correlation 
between  the  models  for  the  0°  soliton  angle. 


Soliton 

Percent  Larger  (AnaIvtical/PE) 

Angle 

Mode  1 

Mode  2 

Mode  4 

Mode  5 

0 

1.55 

5.69 

10.25 

-31.96 

Table  6. 


Maximum  Percent  that  the  Analytical  RME  deviation  from  unity  is 
Greater  than  the  PE  RME  deviation  from  unity 


Soliton  Correlation  Coefficient 

Angle  Mode  1  Mode  2  Mode  4  Mode  5 


0 


0.98806  I  0.88611  |  0.7838  |  0.72149 


Table  7.  Correlation  Coefficients  for  Perturbed/Unperturbed  Ratios  of  Modal 

Energies  at  5  km 


Figure  24.  Eine  Comparison  of  the  Ratio  of  Modal  Energies  at  5  km  (mode  1,  150 

Hz,  0”  soliton  angle) 
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Figure  25.  Line  Comparison  of  the  Ratio  of  Modal  Energies  at  5  km  (mode  5,  150 

Hz,  0”  soliton  angle) 


2,  Comparison  of  the  Analytical  Model  with  the  PE  Model  for  a  Soliton 
Packet  at  Varying  Orientations  to  the  Main  Direction  of  Acoustic 
Propagation 

In  this  section,  the  analytical  model  is  compared  to  the  PE  model  at  varying 
angles  of  the  soliton  packet  by  looking  at  the  surface  plots  of  the  RME.  Eigures  26  and 
27  are  the  surface  plots  of  the  RMEs  for  mode  1,  150  Hz,  at  a  20°  soliton  angle  as 
produced  from  the  PE  model  and  analytical  model,  respectively.  As  in  the  case  of  the  0° 
soliton  angle,  the  20°  soliton  angle  soliton  causes  striation  patterns  in  the  RME  field  after 
the  soliton.  The  patterns  are  believed  to  be  due  to  vertical  mode  coupling  since  the  work 
of  Roush  (2008)  showed  that  horizontal  diffraction  effects  do  not  become  significant  until 
soliton  angles  of  approximately  85°.  Therefore,  for  the  soliton  angles  analyzed  in  this 
paper  the  mode  coupling  effects  will  be  the  result  of  vertical  mode  coupling.  The 
striation  patterns  are  made  up  of  focusing  and  defocusing  regions  in  the  RME  field.  The 
size  of  the  focusing/defocusing  regions  appear  to  be  comparable,  as  well  as  the 
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magnitude  of  the  ratios  in  each  region.  However,  in  the  positive  cross  range  section  of 
Figure  27  there  is  a  focusing  area  that  is  not  present  in  Figure  26.  Again,  the  oscillation 
along  the  focusing/defocusing  regions  is  the  biggest  difference  between  the  two  models. 
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Figure  26.  Ratio  of  Modal  Energies  (mode  1,  150  Hz,  20“  soliton  angle) 
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Figure  27.  Ratio  of  Modal  Energies  (mode  1,  150  Hz,  20“  soliton  angle) 

Next,  the  soliton  paeket  is  rotated  to  50°  to  determine  the  effeet  of  more  grazing 
angles  of  propagation  on  the  eorrelation  of  RME  between  the  models.  Eigures  28  and  29 
are  surface  plots  of  the  RME  for  mode  1,  150  Hz,  and  a  50°  soliton  angle  produced  by  the 
PE  model  and  analytical  model,  respectively.  Both  plots  contain  striation  patterns  of 
comparable  size  and  with  comparable  magnitude  of  RME.  The  oscillations  in  the 
focus ing/defocusing  regions  of  the  PE  plot  are  less  significant  at  this  angle.  This 
observation  casts  some  doubt  on  whether  the  oscillation  is  a  numerical  artifact  due  to  the 
use  of  Cartesian  coordinates  in  the  PE  model.  If  the  use  of  Cartesian  coordinates  was  to 
blame  for  the  oscillations,  then  why  would  the  effect  become  less  significant  with  an 
increased  soliton  angle?  The  2-D  assumption  of  the  theory  is  a  more  likely  source  of  the 
difference,  and  the  oscillations  may  simply  be  an  indication  of  the  horizontal  refraction 
effects.  Erom  Eigures  28  and  29,  we  can  infer  that  the  models  correlate  better  due  to  the 
smaller  effect  of  the  oscillations  in  the  PE  model  with  increasing  soliton  angle. 
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Figure  28.  Ratio  of  Model  Energies  (mode  1,  150  Hz,  50“  soliton  angle) 
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Figure  29.  Ratio  of  Modal  Energies  (mode  1,  150  Hz,  50“  soliton  angle) 
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A  quantitative  look  at  the  RME  will  now  be  done  by  looking  at  them  at  a  fixed 
down  range  value  of  5  km.  Figures  30  and  31  are  the  magnitude  of  the  RME  at  5  km  for 
a  20°  and  50°  soliton  angle,  respectfully.  The  correlation  coefficient  for  the  20°  soliton 
angle  is  0.6678  and  0.4876  for  the  50°  soliton  angle.  This  trend  does  not  support  the 
early  observation  from  the  surface  plots  that  the  plots  correlate  better  for  increasing 
soliton  angles.  In  fact,  values  tabulated  in  Table  8,  which  are  the  correlation  coefficient 
for  varying  soliton  angle  at  each  mode  analyzed,  show  that  there  is  no  clear  trend. 
However,  it  does  show  that  some  conditions  result  in  significant  differences  in  the 
geometric  alignment  of  the  RME. 


Figure  30.  Line  Comparison  of  the  Ratio  of  Modal  Energies  at  5  km  (mode  1,  150 

Hz,  20°  soliton  angle) 
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Figure  31.  Line  Comparison  of  the  Ratio  of  Modal  Energies  at  5  km  (mode  1,  150 

Hz,  50°  soliton  angle) 


Soliton 

Correlation  Coefficient 

Angle 

Mode  1 

Mode  2 

Mode  4 

Mode  5 

0 

0.9881 

0.8861 

0.7838 

0.7215 

10 

0.8796 

0.9230 

0.4275 

0.9666 

20 

0.6678 

0.7524 

0.7295 

0.9214 

30 

0.7897 

0.3533 

0.6867 

0.8111 

40 

0.8318 

0.4524 

0.4485 

0.6482 

50 

0.4876 

0.6988 

0.5664 

0.6794 

60 

0.4887 

0.6765 

0.6342 

0.9375 

70 

0.6025 

0.5699 

0.4101 

0.9366 

80 

0.5497 

0.7476 

0.2828 

0.9058 

Table  8.  Correlation  Coefficients  for  Perturbed/Unperturbed  Ratio  of  Modal 

Energies  at  5  km 


Since  the  correlation  coefficient  only  accounts  for  the  geometric  alignment  of  the 
striation  patterns  and  not  for  any  differences  in  the  magnitude  of  the  RME  within  the 
striation  patterns,  we  have  calculated  the  percent  that  the  RME  deviation  from  unity  is 
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greater  for  the  analytical  model  than  for  the  PE  model  at  the  maximum  difference.  The 
differences  in  the  magnitude  of  the  RMEs  deviation  from  unity  are  tabulated  in  Table  9. 
Again,  these  data  show  that  significant  differences  can  exist  between  the  RME  for  certain 
conditions.  The  analytical  deviation  is  as  much  as  64%  larger  than  the  PE  deviation  for 
mode  5  with  a  40°  soliton  angle.  Eigure  32  is  the  line  comparison  for  mode  5  with  the 
40°  soliton  angle. 


Soliton 

Percent  Larqer  (AnaIvtical/PE) 

Angle 

Mode  1 

Mode  2 

Mode  4 

Mode  5 

0 

1.55 

5.69 

10.25 

31.96 

10 

9.13 

6.60 

4.94 

2.44 

20 

6.93 

7.95 

5.37 

14.80 

30 

7.14 

9.48 

6.07 

5.18 

40 

5.89 

25.29 

10.44 

65.43 

50 

18.49 

26.75 

8.15 

40.25 

60 

26.22 

19.90 

1.84 

28.28 

70 

19.70 

35.42 

6.18 

24.42 

80 

1.80 

35.76 

19.85 

45.18 

Table  9.  Maximum  Percent  that  the  Analytical  RME  Deviation  from  Elnity  is 

Greater  than  the  PE  RME  Deviation  from  Unity 


Eigure  32.  Eine  Comparison  of  Ratio  of  Modal  Energies  at  5  km  (mode  5,  150  Hz, 

40°soliton  angle) 
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3. 


Conclusions 


In  this  section,  we  have  analyzed  the  correlation  between  the  RME  produced  from 
the  results  of  the  analytical  and  PE  models  for  sound  propagation  through  a  soliton 
packet.  Through  those  analyses,  we  have  shown  that  the  analytical  model  does  produces 
results  consistent  with  the  PE  model  over  the  soliton  angles  and  mode  numbers  analyzed. 
Eor  example,  both  models  produce  RME  fields  with  striation  pattern  of  comparable  size, 
shape,  and  magnitude.  However,  some  limitations  have  also  been  discovered.  Eor 
certain  mode  numbers  and  solitons,  the  RME  deviation  from  unity  is  more  than  60% 
greater  for  the  analytical  model  than  for  the  PE  model.  The  correlation  coefficients  for 
the  geometric  alignment  of  the  striation  pattern  also  dropped  to  as  low  as  0.2828.  These 
differences  in  the  results  did  not  have  a  discernible  trend  in  which  to  predict  conditions  in 
which  errors  would  occur.  The  PE  model  also  contained  patterns  in  the  striation  patterns 
that  oscillated  as  a  function  of  range  from  the  source  for  a  given  angle  of  acoustic 
propagation  that  did  not  occur  in  the  analytical  predictions.  These  limitations  make  point 
wise  agreement  of  solutions  unlikely.  However,  the  analytical  model  is  sufficient  at 
predicting  the  same  trends  and  general  behavior  as  produced  by  the  PE  model. 
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V.  SUMMARY  AND  CONCLUSIONS 


In  this  paper,  the  Single  Seattering  Analytieal  Model  is  verified  through  direet 
eomparison  with  the  3-D  implementation  of  a  proven  Parabolie  Equation  model,  the 
3DMMPE.  The  eomparison  in  broken  up  into  three  parts;  aeoustic  propagation  through  a 
single  soliton,  aeoustie  propagation  through  a  soliton  paeket,  and  mode  eoupling 
resonanee  behavior  eaused  by  propagation  through  a  single  soliton.  Eor  eaeh  part,  the 
amount  of  mode  eoupling  predieted  by  eaeh  model  is  used  to  make  the  eomparison.  The 
desired  result  is  not  to  show  that  the  analytie  model  will  produee  point  wise  aeeurate 
results,  but  to  show  that  it  produees  results  eonsistent  with  the  3DMMPE.  Erom  the 
beginning,  the  single  seattering  assumptions  made  in  the  development  of  the  analytieal 
model  were  meant  to  produee  an  estimation  of  the  mode  eoupling  behavior.  The  premise 
of  this  analysis  is  that  the  assumptions  made  produee  results  that  eompare  well  with  the 
3DMMPE. 

Eor  aeoustie  propagation  through  a  single  soliton,  we  show  that  the  analytieal 
model  prediets  equivalent  amounts  of  mode  eoupling.  The  analysis  shows  that  the 
analytieal  model  produees  RME  fields  with  striation  pattern  of  eomparable  size,  shape, 
and  magnitude  to  the  3DMMPE.  The  analysis  also  determines  the  model’s  limitations. 
Eor  eertain  mode  numbers  and  soliton  angles,  the  analytieal  model  predicts  as  much  as 
49%  more  mode  coupling  than  the  PE  model.  The  normalized  correlation  coefficients  for 
the  geometric  alignment  of  the  striation  pattern  can  also  drop  as  low  as  0.7671.  These 
differences  in  the  results  do  not  have  a  discernible  trend  in  which  to  predict  conditions 
where  large  errors  would  occur.  Oscillations  in  the  striation  pattern  present  in  the  RME 
fields  for  the  PE  model  contribute  to  the  differences  seen  in  the  qualitative  analysis  of  the 
RME  fields  and  the  quantitative  comparison  of  magnitude  and  correlation  coefficients. 
These  limitations  make  point  wise  agreement  of  solutions  unlikely.  However,  despite 
this  limitation,  the  analytical  model  is  sufficient  at  predicting  the  same  trends  and  general 
behavior  as  produced  by  the  3DMMPE  model. 

Two  more  solitons  were  added  to  the  acoustic  field  to  determine  its  affect  on  the 
analytical  model’s  ability  to  produce  consistent  predictions  for  acoustic  propagation 
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through  a  soliton  packet.  In  the  analysis,  we  show  that  the  analytieal  model  still  prediets 
equivalent  amounts  of  mode  eoupling.  While  the  predietions  are  not  as  elose  as  for  the 
single  soliton,  the  model  still  produces  RME  fields  with  striation  pattern  of  comparable 
size,  shape,  and  magnitude  to  the  3DMMPE  model.  Again,  limitations  of  the  analytieal 
model  were  identified  by  the  analysis.  Eor  eertain  mode  numbers  and  solitons,  the 
amount  of  mode  eoupling  predieted  by  the  analytieal  model  was  as  mueh  as  60%  more 
than  what  was  predieted  by  the  PE  model.  The  eorrelation  eoeffieients  for  the  geometric 
alignment  of  the  striation  pattern  also  dropped  to  as  low  as  0.2828.  These  differences  in 
the  results  did  not  have  a  diseernible  trend  in  whieh  to  prediet  eonditions  in  whieh  errors 
would  oecur.  The  oseillations  in  the  striation  patterns  present  in  the  RME  field  for  the  PE 
model  eontinued  to  aecount  for  significant  differences  between  the  models  shown  by  the 
analysis.  Again,  these  limitations  make  point  wise  agreement  of  solutions  unlikely,  but 
the  analysis  shows  that  the  analytieal  model  is  still  suffieient  at  predicting  the  same  trends 
and  general  behavior  as  produeed  by  the  PE  model. 

The  mode  eoupling  resonanee  behavior  predieted  by  the  analytieal  model  is  also 
eonsistent  with  the  observed  resonanees  in  the  PE  model  predictions.  The  analysis  shows 
that  the  analytieal  model  ean  prediet  which  soliton  width  the  PE  model  will  produee  a 
resonanee  for  within  17%  error  for  mode  1.  Eor  mode  2,  the  analytieal  model  is  less 
aeeurate  at  predieting  the  resonanees  of  the  PE  model.  The  errors  in  the  resonances  are 
as  high  as  35%,  with  one  outlier  as  high  as  70%.  The  inconsisteney  in  the  predietion  of 
the  resonances  for  the  different  modes  makes  it  elear  that  more  eomparisons  are  needed 
to  determine  just  how  well  the  analytical  model  can  predict  the  resonances.  Further  analysis 
should  be  conducted  to  include  more  modes  and  the  reproduction  of  Figure  17  using  the 
correct  initial  conditions. 

Further  researeh  should  also  be  condueted  to  determine  the  souree  of  the 
oseillations  in  the  striation  patterns  present  in  the  RME  fields  for  the  3DMMPE  model. 
The  ehange  in  the  oseillation  with  mode  number  and  soliton  angle  suggest  that  they  are  a 
physieal  phenomena  not  aeeounted  for  by  the  2-D,  single  seattering  approximation.  One 
possible  explanation  eould  be  that  they  are  eaused  by  3-D  horizontal  refraetion  effects.  A 
comparison  of  the  3-D  MMPE  and  the  2-D  MMPE  eould  resolve  the  oseillation  to  3-D 
effeets. 
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